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Abstract 

A field theoretic renormalization group method is presented which is capable of deal- 
ing with crossover problems associated with a change in the upper critical dimension. 
The method leads to flow functions for the parameters and coupling constants of 
the model which depend on the set of parameters which characterize the fixed point 
landscape of the underlying problem. Similar to Nelson's trajectory integral method 
any vertex function can be expressed as a line integral along a renormalization group 
trajectory, which in the field theoretic formulation are given by the characteristics 
of the corresponding Callan-Symanzik equation. The field theoretic renormalization 
automatically leads to a separation of the regular and singular parts of all crossover 
scaling functions. 

The method is exemplified for the crossover problem in magnetic phase transi- 
tions, percolation problems and quantum phase transitions. The broad applicability 
of the method is emphasized. 



1 Introduction 



Systems in the vicinity of a critical point show scale invariance [1]. Quantities 
like the susceptibility obey a homogeneity property, which is characterized by 
critical exponents and scaling functions. Within the framewrk of renormaliza- 
tion group theory it can be shown that this scaling property is closely related 
to the existence of a fixed point of the renormalization group transformation. 
In the simplest case the critical behavior of a particular problem is dominated 
by just one fixed point. Quite frequently, however, one encounters physical 
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systems which can exhibit different scaling behavior in different asymptotic 
regimes. In those cases it becomes important to study the crossover from one 
type of critical behavior to another. Examples for such crossover problems 
are bicritical and tricritical points (in general multicritical points), magnetic 
systems with anisotropic interactions such as anisotropies in the exchange in- 
teraction as well as dipolar interaction, finite size problems, and in general 
any problem which is characterized by some length scales in addition to the 
correlation length £. 

The theoretical investigations of those crossover phenomena started with a 
scaling theory by Riedel and Wegner [2], and the extended scaling hypothesis of 
Pfeuty, Fisher and Jasnow [3] . Nelson has presented the first renormalization- 
group formalism for calculating the crossover scaling functions associated with 
such situations, the so called trajectory integral method [4]. In this paper we 
will present a field theoretic version of Nelson's trajectory integral method 
and exemplify the method by several problems from various fields of critical 
phenomena. The method is based on the work by Amit and Goldschmidt [5] 
who proposed a "generalized minimal subtraction method", which has sub- 
sequently been applied to various crossover problems [7,8,11]. Recently the 
method has been extended to deal with situations, where the crossover is 
accompanied by a change in the upper critical dimension [9,10]. 

In many crossover problems one encounters a situation where besides the cor- 
relation length £ there is a second length scale, characterized by the anisotropy 
scale g. Out of many examples we would like to discuss three quite character- 
istic cases. 



(1) The critical behavior of a uniaxial dipolar ferromagnet is described by the 
effective Ginzburg-Landau- Wilson (GLW) free energy functional [9] 
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where g measures the strength of the dipolar interaction relative to the ex- 
change interaction, and the wave vector k has been decomposed into a com- 
ponent q along the uniaxial direction and the remaining components p. 

(2) The crossover from isotropic to directed percolation can be described by 
a type of Reggeon field theory with the action [12,13] 
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(3) Multicritical behavior near the superfluid to Mott insulator quantum 
phase transition in Bose lattice models are described by the effective ac- 
tion— citeFGW89 

5 = 1 J (r + k 2 + 2giu + u 2 ) | 0(k, u) | 2 +^ / |0(x,t)| 4 , (3) 

where the field serves as an order parameter for superfluidity and the param- 
eter g (which vanishes at the multicritical point) is a measure of the absence 
of particle-hole symmetry. 

All these phenomena have in common the characteristic that the anisotropy 
scale appears in the harmonic part of the effective functional of the model, and 
that asymptotically there is a change in the upper critical dimension (usually 
a reduction). Hence there are not only different critical points complicating 
the renormalization group analysis, but in addition those fixed points are 
characterized by different upper critical dimensions. The latter invalidates 
conventional methods such as ^-expansion since there is no fixed upper critical 
dimension to expand around. On the contrary the crossover may be even 
viewed as a continuous variation of the upper critical dimension. 

The paper is organized as follows. In the following section we introduce an ex- 
tended minimal subtraction procedure which is capable of describing crossover 
phenomena from one fixed point to another accompanied by a change in the 
upper critical dimension. In section 3 we exemplify the method for the three 
models mentioned above. Finally, we summarize our results in section 4 and 
give an outlook concerning the broad applicability of the method. 



2 The method of generalized minimal subtraction 

Crossover phenomena are characterized by the presence of at least one addi- 
tional relevant length scale g besides the correlation length £ oc r~ v '. Frequently 
it is given by an anisotropy or "mass" parameter describing the variation from 
one scaling region to another, or in renormalization group language, from a 
primary fixed point to a secondary fixed pointQ. We have listed some exam- 
ples in the introduction. In order to describe this crossover from one critical 
behavior to another Riedel and Wegner [2] have formulated a scaling theory, 
and later Pfeuty, Fisher and Jasnow [3] proposed an extended scaling hypoth- 
esis for crossover systems. Let r be the reduced temperature and g be the 

2 In general the scenario may be even more complicated. Depending on the number 
of length scales present in the problem there may by more than two fixed points 
with a cascade of various crossover phenomena. 
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scaling field characterizing the additional length scale in the problem, then 
the crossover scaling Ansatz, e.g. for the susceptibility, reads 



X(r,g)=r-if(g/r*) 



(4) 



near the multicritical point (r, g) = (0, 0). Here 7 is the susceptibility exponent 
at the primary fixed point, while <fi is the crossover exponent. In order to 
describe the change of the critical indices induced by the crossover from the 
primary to the secondary fixed point, one assumes that the crossover scaling 
function f(x) becomes singular at x c 



as r — > r c (g) approaches criticality for fixed g. As is obvious from the above 
description of the phenomenological crossover scaling theory, it has the disad- 
vantage that the crossover in the critical indices is incorporated in the crossover 
scaling function in a rather singular way. It would be much more convenient 
to have a formulation of the crossover problem, where the crossover of the 
exponents is incorporated in the power law prefactor r~ 7 by 7 becoming an 
effective exponent depending on the crossover scaling variable x — g/r^. This 
would leave us with a regular scaling function f'(x). As will become clear 
later the field theoretic method proposed in this section, will satisfy the latter 
requirement. 

The first quantitative theory for crossover scaling functions based on the renor- 
malization group theory was given by Nelson [4]. He introduced a formalism 
which expresses the free energy as a trajectory integral along the renormal- 
ization group flow lines. The method proposed below will be similar in spirit 
to this trajectory integral method, but will have considerable technical ad- 
vantages due to the explicit separation of singular and regular terms in the 
crossover scaling function. 

Now we describe a field theoretic method for the analysis of crossover phenom- 
ena. The technical difficulty one encounters here is that both the ultraviolet 
(UV) and infrared (IR) singularities will differ in the two distinct scaling 
regimes. Early work on formulating crossover phenomena in terms of field 
theoretic methods, tried to parallel the formulation using phenomenological 
scaling theories. In this "traditional" approach to crossover problems, one 
would compute the critical exponents at one of the stable fixed points; all the 
crossover features would then be contained in the accompanying scaling func- 
tion as corrections to this scaling behavior. However, in general a calculation 
to high order in the perturbation expansion would be required in order to 
achieve a satisfactory description of the entire crossover region. Of course, us- 
ing an (e = d c — d) expansion with respect to either of the fixed points renders 
the other one completely inaccessible, if their upper critical dimensions do 
not coincide. Therefore Amit and Goldschmidt's idea [5,6] to incorporate the 
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crossover features already in the exponent functions has proven much more 
successful than treating the problem on the basis of scaling functions. The 
essential prescription one has to bear in mind is that the renormalization con- 
stants are not solely functions of the anharmonic coupling, but necessarily also 
of the additional "mass" or anisotropy parameter g Q describing the interplay 
between the two different scaling regimes [5]. For a consistent treatment of the 
entire crossover region, one thus has to assure that the UV singularities are 
absorbed into the Z factors for any arbitrary value of g , including g — > oo. 
This is not a trivial prescription, as usual the l/(d c — d) poles will be altered 
in the different scaling regimes. For the situation that we have in mind, even 
the value of the upper critical dimension is bound to change as the crossover 
takes place, in contrast to previously studied cases [7,8,11]. However, we shall 
demonstrate that with the above stated so-called "generalized subtraction 
scheme" this change in the upper critical dimension may be incorporated into 
the usual formalism without any drastic changes, provided one refrains from 
any e expansion about d c . The perturbation series is then an expansion with 
respect to the effective coupling v to be introduced later, which is not an a- 
priori small parameter, and the perturbation expansion is uncontrolled in this 
sense. If higher orders of the perturbation expansion were known, one could 
substantially refine the theory by a Borel resummation procedure. For a more 
detailed discussion of the question in which cases one may dispense with a 
(d c — d) expansion, we refer to work of Schloms and Dohm [15]. 

We remark that the somewhat misleading term "generalized minimal subtrac- 
tion scheme" stems from the fact that in the framework of an e expansion this 
corresponds to adding logarithms of go, which are finite in the limit e — > 0, 
to the Z factors. This is the procedure proposed in the original work by Amit 
and Goldschmidt [5]. But it works only if the upper critical dimension of the 
primary and secondary fixed point are equal. Otherwise the method by Amit 
and Goldschmidt has to be extended. There are several ways to implement the 
requirement that the vertex functions become finite in both limits go — > and 
go — > oo after reparameterizing the theory using renormalization factors. We 
will explicitly use two of them. Firstly, for the uniaxial dipolar ferromagnet, 
discussed in section 3.1, we will use the prescription of Amit and Goldschmidt, 
extended in the following way. Upon performing an e-expansion around the 
upper critical dimension of the primary fixed point, one gets to one loop order 
renormalization factors of the structure Z = 1 + - + bin go. We reexponentiate 
those logarithms of g such as to map the result one would obtain by doing an 
expansion of the model for g — > oo around the upper critical dimension of the 
secondary fixed point. This prescription allows for an analytic treatment of the 
flow equations. Secondly, we will determine the renormalization factors by re- 
quiring that it contains all the singularities for both limits go = and go = oo, 
but without performing any expansion around any of the critical dimensions. 
A (minor) price to be paid is that for the flow equations and related quanti- 
ties only numerical solutions are accessible, and merely the limiting cases of 
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go = and go — > oo, respectively, allow for an analytical investigation. The 
latter method is quite close to the method of normalization conditions [6], but 
is minimal in that it contains just the poles in ^ for both limits go = and 
g — > oo. Here e refers to the distance to the upper critical dimension of the 
primary and secondary fixed point, respectively. 



2.1 Flow equations 



Although the following discussion of the flow equations can be presented more 
generally, it is convenient to have the specific set of effective free energy func- 
tionals and actions in mind, which we have presented in the introduction. 
Therefore, we assume that all of Wilson's flow functions characterizing the 
crossover in the fixed points and the corresponding critical exponents depend 
only on two parameters, ( = ((u,g). 

The renormalization group equations relate the values of an arbitrary vertex 
function at one point in parameter space to a transformed point. Conceptually 
the renormalization group equations are obtained by observing that all bare 
vertex functions T are independent of the arbitrary renormalization scale \x 

/x|-| o r (a,^, M ;q) = 0, (6) 

where a is a set of parameters characterizing the effective free energy functional 
under consideration, and | indicates that all derivatives have to be taken at 
fixed values of the bare parameters. The quantity g is the second length scale 
of the problem besides the correlation length £, and u is the coupling constant 
of the nonlinearity. Wilson's flow functions are defined by 

d 

Cx = ii-^\ ln(x/xo) , (7) 

where x stands for the parameters a, the anisotropy scale g and the coupling 
constant u. Furthermore we introduce the (^-function for the renormalization 
of the vertex function Y 

d 

Cr = /^lo ln ^r, (8) 
where T = Z-pT. Then one derives the following partial differential equation 



r(/x, a, g, u; q) = (9) 
for the renormalized vertex function. As a consequence of the generalized 
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renormalization scheme, all the flow functions are functions of the coupling 
constant u as well as the anisotropy scale g. This leads to renormalization 
group trajectories in a two-dimensional parameter space, the (g, w)-plane. 
Within this plane there are usually several fixed points, where two of them 
are of particular importance. These are the primary fixed point P located at 
g = and the secondary fixed point S located at g — oo. The primary fixed 
point turns out to be infrared stable for g = only, whereas all the renor- 
malization group flow tends to the secondary fixed point S if the anisotropy 
scale g is nonzero. The renormalization group flow within Wilson's elimina- 
tion procedure is generated by the introduction of a spatial rescaling factor 
b. In the present field theoretic framework, a renormalization group trajec- 
tory is given in terms of the characteristics of the Callan-Symanzik equation, 
Eq.(9). The characteristics are defined as the solution of the following first 
order differential equations 

^ = CW) > (io) 

with the initial conditions x(l) = x, namely x(£) = xexp ( x (£')d£' /£' . If we 
denote the engineering dimension of the vertex function T by d-p and define 
the dimensionless (renormalized) vertex function T by 

r( / i,{x};q,^) = /i dr f({x};q/ yU ), (11) 

the solution of Eq. (9) reads 

r(/i,{x};q) = (/tf)*exp j/ ^ylf^^uMjij • (12) 

Here we have introduced an effective anharmonic coupling v, which is some 
combination of the nonlinearity u and the anisotropy scale g, where the de- 
tailed form depends on the particular problem one is considering. The intro- 
duction of this effective coupling constant is necessitated by the fact that one 
should work with coupling constants whose fixed point values are finite in 
both limits, g — > and g — > oo. The representation (12) of the vertex func- 
tion in terms of the characteristics of the parameters and coupling constants 
constitutes an explicit expression for the desired crossover scaling function^]. 

Compared to the trajectory integral method of Nelson [4] the above formalism 
has the advantage that the singular and the regular part of the crossover 
scaling functions are explicitly separated. This is simply a consequence of the 
field theoretic renormalization procedure. 



3 Here we do not consider the case of composite-operator vertex functions, which 
have to be renormalized additively. For an extension of the method to include these 
additively renormalized quantities we refer the reader in particular to Ref. [8] 
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2.2 Flow diagram, fixed points and exponents 



In this subsection we give a qualitative discussion of the renormalization group 
flow in the (g, i>)-plane, where g = g/(l + g) and v are appropriately defined 
effective coupling constants, which are finite in the isotropic (g = 0) and 
anisotropic (g = oo) limit. Usually the topology of the flow diagram is de- 
termined by the presence of four fixed points, the primary and the secondary 
fixed points, and two Gaussian fixed points with v — at g — and g = oo. 
The primary fixed point is infrared stable only if the anisotropy parameter 
equals zero. Otherwise the only infrared stable fixed point is the secondary 
fixed point. All other fixed points are unstable. The flow diagram is separated 
into two regions by a separatrix, which is given by the trajectory originating 
from the primary fixed point and terminating in the secondary fixed point. 
This renormalization group trajectory describes the universal crossover from 
primary to secondary critical behavior. 

The flow diagram can be investigated quantitatively upon solving the flow 
equation of the running coupling v(£) 

Mt) = PM , (13) 



d£ 



v. In the flow equations above, 

o 



where the (3 function is defined by (3 V = fj,-^- 

the parameter £ may be considered as describing the effect of a scaling trans- 
formation upon the system. Roughly speaking In £ corresponds to the Wilson 
rescaling factor b. Obviously, the theory becomes scale-invariant when a fixed- 
point v*, to be obtained as a zero of the (5 function, /3 v (v*) = 0, is approached. 
The properties of the vertex function in the vicinity of the fixed point will 
yield the correct asymptotic behavior, if the latter is infrared-stable. Upon 
defining the fixed point values of the flow functions by Q = ( a (v = v*), also 
called the anomalous dimensions of the parameters a, we find 

r(/i, {x};q) = /Af^f </*V; ^ . (14) 

Using the matching condition £ = q/fi, one arrives at the following scaling 
form 

r(jM, r, g, v*; q) oc q 2 ^t u (jjtffiQZ> v*; , (15) 

where we have specialized to the simple case that the set of parameters a, 
reduces to just one parameter, namely the reduced temperature r. We defined 
two independent critical exponents according to r\ = —Q and v = — ^. The 
crossover exponent <p can be read off from the homogeneity relation for the 
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vertex function upon using the matching condition r(£) oc rft* = 1 

r>, r, g, v; q) oc rT (l, g/r*, v*, ; V") , (16) 

where 7 = z/(2 — 77) and = Q/C- Note that all of the above exponents 
correspond to the exponents at the secondary, and asymptotically stable, fixed 
point. If one could neglect the variation of the scaling functions on the flow 
of the parameters on the right hand side of Eq.(12), the effective exponents 
of the problems would simply be given by equations like rj = where the 
fixed point values of the Wilson flow functions are replaced by their values 
along the renormalization group trajectory. Such an approximation is quite 
frequently referred to as a "renormalized mean field" approximation. 



3 Model systems with dimensional crossover 

3.1 Uniaxial dipolar j err omagnet 

The influence of the dipolar interactions on the critical behavior of magnetic 
systems with isotropic and uniaxial exchange interaction is quite different, 
especially if one takes into account the dipole-dipole interaction present in any 
real ferromagnetic material. For isotropic ferromagnets the dipolar interaction 
leads only to a slight modification of the critical exponents. Uniaxial dipolar 
ferromagnets, however, show classical behavior with logarithmic corrections 
in three dimensions [16-18], which in the language of renormalization group 
theory means that the upper critical dimension for the asymptotic critical 
behavior is shifted to d c — 3. 

The existence of logarithmic corrections was verified experimentally for a num- 
ber of uniaxial ferromagnetic substances [19,20]. However, these experiments 
were performed in regions of reduced temperature, where departures from the 
asymptotic behavior are expected and are indeed observed. In particular, one 
finds a maximum in the effective exponent of the susceptibility [20]. Here we 
study this crossover by a the extended minimal subtraction described in the 
preceding section. This method allows the calculation of the complete flow 
of the coupling constants and parameters from the Ising fixed point to the 
uniaxial dipolar fixed point. 

The combined effect of the uniaxial anisotropy of the exchange interaction 
and the anisotropy of the long-range dipolar interaction leads to the following 
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effective free energy functional [16-18,21] 



H 



k L 



r + k 2 + g 2 — 



0«(k)0-(-k) + ^/[0«(x)0«(x)] 2 . (17) 



This problem is characterized by two length scales. Firstly, there is the cor- 
relation length £, which diverges £ oc r~ u as the critical temperature T c is 
approached. Here r is the renormalized value of the bare reduced temperature 
variable r = (T — T c °)/T c °. Secondly, there is a length scale introduced by the 
presence of the dipolar interaction, given by the so called dipolar wave vector 
Id oc g$. 

The renormalized parameters, coupling constants and fields are defined by 



u 



fi d 4 Z u S d u , r = Z^ZrTon 2 , 9 = Z^Z g g^\i 1 , and ^g(k) = Zj a (k), 
where the factor S d = 2/(4vr) d / 2 r(rf/2) is introduced for convenience. One 
finds to one loop order 



Z, k = 1 



u — (1+9) 



Z r = 1 



n + 2. 



(18) 



and Z^ — Zg — l. The susceptibility can be analyzed in terms of the effective 
exponent 7 e fj = <ilnx _1 (r, g, u)/d\nr. The one loop result is 

d\n X ^(l,g(i),u(e)) 



7 eff = -2/CrW),«W) + 



dlnr 



(19) 



with the matching condition r(£)/n 2 £ 2 = 1. The ^-function for the mass, is 
C r (£) = -2 + {n + 2)w(f)(l + #(f))-76, and for d = 3 the second term in 
Eq.(19) reads 



d\n X ~\l,g(£),u(£)) 



dlnr 



1/2 



X 



l+(l+g- 1 (tj)~ 1 -2g(e) \n(l+g-\£)) 



(20) 



The solution of the flow equation, £du(£)/d£ = (5{£) with (5 = /id^u\ , for e — 1 
is 



u{£) = u H g{£) 



+ In 



1+9 



(21) 



where w# = 6e/ (n + 8) is the Heisenberg fixed point value. In the asymptotic 
limit the susceptibility reduces to x l { r ) — r I mr |~( n+2 )/( n+8 ) in agreement 
with Refs. [17,18] In Fig.l the experimental data from Ref. [20] for 7 e g- are 
compared with the theory, where the non universal parameters in Eqs.(19)- 
(21) have been chosen as g — 1 and u = u H = 1/3 (n = 1). In the range r < 1 
of the reduced temperature there is excellent agreement between theory and 
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Fig. 1. Effective exponent of the susceptibility for g = 1.0 and u = u 7 versus 
ln(r). The momentum scale is chosen as p = 1. Experimental data on LiTbF^ from 
Ref.[20]. 

experiment. Especially, it is found that the observed crossover corresponds to 
the flow from the Ising fixed point to the uniaxial dipolar fixed point. For r > 1 
the data tend to the mean field value %g = 1 corresponding to a crossover 
to the Gaussian fixed point. This can be described only qualitatively within a 
critical theory by the limit £ — + oo. 



3.2 Crossover from isotropic to directed percolation 

Another example, where the presence of a second length scale leads to a 
crossover in the critical behavior, is the physics of biased percolation. The 
patterns generated by these percolation models incorporate self-similar as well 
as self-affine clusters. 

In ordinary percolation, sites or bonds are filled at random with probabil- 
ity p [22]. The percolation process then proceeds along paths connecting oc- 
cupied nearest neighbors. The clusters formed by nearest-neighbor links are 
self-similar, i.e., they display isotropic scaling. In directed percolation [23] the 
links between nearest neighbors have a bias in one preferred direction t, such 
that the percolation process advances along this direction only. The size of 
the clusters in the preferred direction is characterized by a length scale dif- 
ferent from that in the perpendicular direction. If percolation in the positive 
t direction is merely favored with a certain probability with respect to the 
negative t direction, but propagation "backward" in time is still admitted, the 
situation will be more complicated. If the "anisotropy" g is low, one expects 
almost isotropic scaling behavior in a large region of the phase diagram. How- 
ever, when the critical region near the percolation threshold p c is approached, 
self-affine scaling will become apparent. 
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In order to give a quantitative description of the crossover from isotropic 
to directed percolation we investigate the pair correlation function G(r 2 ,ri), 
which measures the probability that sites r 2 = (x.2,t 2 ) and r 1 = (xi,£i) 
belong to the same cluster. Following the work of Cardy and Sugar [12], 
and Benzoni and Cardy [13], the pair correlation function (the superscript 
"0" denotes unrenormalized quantities) is found to be a sum G°(r2;ri) = 
Emn=i[( _ *) m+ "~ 2 / m ' n ']Cmn( r 2; r i)! where the correlation functions are 
obtained from the effective action 



J =J d D xjdt^(p (x,t) 



uo 
2 



r -V 2 -\d 2 + 2 -^d t 
Cq c 

4> (X, t) 2 O (x, t) - O (X, t)0 O (x, t) 



0o(x,t) 

(221 



with r oc p — p c . The effective action is renormalized upon introducing renor- 
malized fields <p = Z^ 2 4>o an d = Z^ 2 (f) , dimensionless renormalized pa- 
rameters r = Z^Zrfo - r 0c )fi' 2 , c 2 = Z^Z' 1 ^, g = Z^ 1/2 Z~ 1/2 Z^^' 1 , 
and a renormalized coupling constant u = Z^ 3 ^ 2 Z u UqB 1 J 2 ] p( d ~ 6 )/ 2 . B d = T(4 — 
d/2)/(47r) d//2 is a geometric factor, and r 0c denotes the fluctuation-induced 
shift of the percolation threshold. The five independent multiplicative renor- 
malization constants are then uniquely determined by extracting the ultravi- 
olet singularities of the two- and three-point vertex functions at finite mass 
r = fi 2 . In the course of renormalizing the model it is convenient to introduce 
an effective anharmonic couplingQf = u 2 clf 7 (g). 

From an analysis of the renormalization group equation (where for notational 
convenience we have taken \x = 1 in the following) 

r u ( W , q , u) =/» e Jr^/'f 11 (r M , B «),S ^.^) . (23 ) 

one finds four independent critical exponents. At a fixed point v* using the 
matching condition I = q, one arrives at the following general self-affine scaling 
form 

r u (r, c, g, u, q, u) oc g—f n (-^, v*, 1, , (24) 

with the critical exponents r/± = —(1, v±_ = z = 1 + Q — Q, and 

z A = — Q. Here, r]± and v± correspond to the two independent indices familiar 
from the theory of static critical phenomena. The exponent z was introduced 
in analogy to a dynamical critical exponent, and in our case is related to 



4 Here I^ n {g) = Jq x m l 2 1 (l + xg 2 )^ d m n ^ 2 dx (with odd integers m and n) de- 
notes a certain type of integrals, with J^ n (0) = 2/m, and lim s _> 00 [g m /^ m (<7)] = 
T(m/2)T((n-d)/2)/T{(m + n-d)/2). 



12 



0.0 0.2 0.4 0.6 0.8 1.0 

gd)/[l+g0)] 



Fig. 2. Renormalization-group trajectories in the (v,g) parameter space for the 
crossover from isotropic to directed percolation (to one-loop order) in 6? = 3 dimen- 
sions. 

the anisotropic scaling behavior. Finally, A is a positive crossover exponent 
describing the transition from isotropic to directed percolation. It stems from 
the fact that there appear two different scaling variables for the "frequency" 
uj in Eq. (24). In the asymptotic limit of directed percolation, g — > oo, the 
second scaling variable disappears, and the scaling behavior is described by 
the three exponents rj±, u±, and z. 

Similarly, with the choices £ = (gu / c) 1 ^ 1+ '* c ^ and £ = r -1 /^ the longitu- 
dinal exponents 2 — r}\\ = (2 — r}±)/z and u\\ = zu± and the ("susceptibil- 
ity") exponent 7 = v±(2 — r)±) = v\\(2 — r)\\), respectively. In the special case 
of isotropic percolation (go = 0), the scaling relations simplify considerably, 
describing self-similar scaling with only two independent critical exponents 
V = -Q and v = -1/C;. 

We have seen that both the self-similar and the self-affme scaling behavior 
are within the scope of the present theory, at least for dimensions d < 5; for 
5 < d < 6 the model is not renormalizable in the directed limit, and sim- 
ply characterized by the exponents corresponding to the Gaussian fixed point 
Vq D , with logarithmic corrections for d = 5. This again emphasizes the fact 
that no expansion with respect to a fixed upper critical dimension can be ap- 
plied consistently. In order to study the entire crossover region between these 
asymptotic regimes, the coupled set of flow equations have to be solved nu- 
merically. Fig. 2 shows the flow diagram (at d — 3 dimensions) for the effective 
couplings v(£) and g(£), whose topology is determined by the four fixed points. 
The separatrix connecting v{ with the infrared-stable fixed point v^, describes 
the universal crossover features from isotropic to directed percolation. 

The interchange from self-similar to self-affine scaling is most conveniently 
described by introducing effective exponents for the pair-correlation function. 
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Using the zero-loop result for the two-point vertex function 



ru({iUw) = ^ ww 

we specialize to r = and q = and define 



2 c(£) 2 £c(£) 



(25) 



din J | rn(o,o,w) 



2 



2 dl^ ■ < 26 > 

Using |cj 2 /£ 2 c(£) 2 + 2^(£)/£c(£)| 2 = 1 yields 2-^, cfr (£) = [2+C^(£)]rfln£/dlnc;. 
Similarly, in the case r = u = one finds 

rfln 1^(0,^0) 

2 - ?7±efT(?) = jj , (27) 

amq 

which reduces to 2 — i]j_ c g = 2 + (,$(&), if {q/£) 2 = 1 is inserted. Finally, 
considering q = and uj = one introduces 

rflnr n (r,0,0) 

7e S W = ^ , (28) 

and choosing the matching condition r(£) = 1 we find 7 e ff (£) = — [2+( ( p(£)]/( r (£). 

The flow of the effective exponents r]\\ e s(£), f]± c s{£), and / j e s{£) in d = 3 
dimensions is depicted in Fig. 3, with the initial value for the coupling v(l) — 
Vj of the isotropic scaling fixed point. The dependence on the anisotropy scale 
g was eliminated by plotting versus the scaling variable lng(£); the graphs 
corresponding to different initial values g(l) then all collapse onto one master 
curve. The most important conclusion to be drawn from Fig. 3 is that the 
anisotropy scale, at which the crossover occurs, considerably differs for the 
effective exponents defined above. r}\\ e g starts to cross over from the isotropic to 
the directed fixed point value already at In g(4ross) ~ —0.8, whereas 7 c g shows 
this crossover at lng(£ cross ) ss —0.2, and ??_i_ c ff only at lng(£ cross ) +0.8. Note 
that the sizeable change of 7711 e s is already apparent at mean-field level, where 
it acquires the values and 1 in the isotropic and directed limit, respectively. 
However, a crossover of the exponents ?7_i_ e ff and j e s requires the ( functions 
at least on the one- loop level. 



3.3 Multicritical behavior in lattice boson models 



The critical behavior at the zero temperature Mott insulator to superfluid (S- 
I) transition of lattice boson models can be analyzed in terms of the effective 
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Fig. 3. Effective exponents r]^ eS (£) (solid), g± e s(^) (dashed), and 7 e ff(^) (dotted) for 
the connectivity as functions of the running anisotropy parameter g(£) for a fixed 
initial value of the coupling constant v(l) = vf. 



action [14] 



r + k 2 + uj 2 + 2ig u \ <f) (k, uj) \ 2 + 



u 



0o(x,t) 



(29) 



where the field serves as the order parameter for superfluidity. The parame- 
ters r and u can be calculated from the underlying microscopic lattice mod- 
els. One finds that go(fi, J) = <9r (/i, J)/d/i, where \i and J are the chemical 
potential and the hopping matrix element of the lattice model, respectively. 
Since in mean field theory the phase boundary is given by ro(/x, J) = 0, the 
point go — constitutes a multicritical point for S-I transition. 

Before starting the renormalization group calculation of the crossover scal- 
ing functions, we give a brief review of the phenomenological scaling the- 
ory [14]. The quantum phase transition at T = is characterized by both a 
diverging length scale £ oc r~ v and a diverging time scale fi" 1 oc r~ vz . Near 
the multicritical point (r, g) = (0, 0) the singular part of the free energy shows 
crossover scaling behavior 



f s Kr^f s (g/r*), 



(30) 



with being the crossover scaling exponent. As discussed in section 2 the 
crossover scaling function f s (x) has to become singular at some critical value x c 
in order to describe the crossover from primary to secondary critical exponents. 



Here we use the field theoretic method described in section 2 to determine 

1/2 

the effective exponents. We introduce renormalization factors by = 0o, 
and dimensionless renormalized parameters r = Z^ x Z r {ro — r-Q C )n~ 2 , g = 

Z^ 2 Zgg^yT 1 , and a renormalized coupling constant u = Z^ 2 Z u u B d fi d ~ 3 . 
The crossover induced by the term linear in the frequency u leads to a reduc- 
tion of the upper critical dimension from d c — 3 to d c — 2. From a one-loop 
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calculation we find [24] for the renormalization factors — 1, Z g — 1, and 



i 

Z r = 1 



^ J dxx m ' 2 ~\x + g^-m-n)/2 ^ ^ 
o 

Z u = 1 - - , (32) 

where A = \lf 2 {.g) and B — (1 + g 2 ) 6 ^ 2 with 5 = 3 — c? being the distance 
to the upper critical dimension of the primary fixed point. One should note 
that Z g = 1 is valid to all orders in perturbation theory since there is no 
singular contribution proportional to uj. The one-loop theory may be further 
improved by taking into account results which become exact in the limit g — > 
oo. In this limit the special ("causal") structure of the propagator implies 
that the corrections to the bare two-point vertex function vanish to all orders 
in perturbation theory [25,14], i.e., Z r (g — > oo) = 1 exactly. Furthermore, 
there are only "ladder" -diagrams [25] contributing to the renormalization of 
u, which form a geometric series with the result 

e M / oo \ 

w = TT^ ' (33) 

with A = ^BdH~ e A. Now, in order to improve the one-loop result we suggest 
to resum the ladder diagrams corresponding to the one-loop term Auo, while 
still using the one-loop result for all the other terms. This procedure amounts 
to 

ufi ' = TTJn, ~ Bu ° ' <34) 

where B = -B d [i~ e B. Equivalently to one loop order we could also use ufi £ = 
u /(l + (A + B)uo). 

Upon defining Wilson's flow functions as in section 3.2 one can study the be- 
havior of the vertex function r n under renormalization group transformations 

r n (/i, r, g, u; q, oo) = ft 2 f n (r(£), g(£), v(£); ± , (35) 

where it is convenient to introduce the effective coupling constant v = ulf 2 (g). 
At the secondary fixed point v* one finds using the matching condition r(£) = 1 

Tn(/i, r, g, u; q, w) oc £ 2u T n (l, g£, v*; q£, uj£) , (36) 

where £ oc r~ u . Since Z g — 1 to all orders in perturbation theory, one finds 
that the crossover exponent <fi is exactly given by the correlation length expo- 
nent. One can define an effective exponent v{f) = —l/( r (£), the flow of which 
is obtained from a solution of the flow equations [24] . An effective dynamic ex- 
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ponent may be obtained by analyzing the pole T of the dynamic susceptibility 
(we have used the 0-loop result for the scaling function) 



r oc frt g(£) - J, 



r{£) + g 2 (i) + q 2 /(fi£) 2 



(37) 



the scaling behavior of which may be analyzed in various ways [24]. At the 
critical point r = with q/fj,£— 1 we get 



describing the trivial crossover from z = ltoz = 2asq^0. Note, that there 
is, however, nontrivial crossover behavior for the exponent of the correlation 
length from the Ising fixed point value v\ to the Gaussian fixed point value 



4 Summary and outlook 

In this paper a field theoretic renormalization group technique using normal- 
ization conditions was presented, which allows for the description of crossover 
phenomena, where the primary and secondary fixed point have different upper 
critical dimension. This method can be thought of as the field theoretic version 
of Nelson's trajectory integral method. It has the advantage of being capable 
of describing quite complicated crossover scenarios. Especially, using the field 
theoretic renormalization method one can derive an explicit expression for the 
crossover scaling function, where the regular and singular parts are already 
separated. 

We have employed this method for three particular interesting crossover phe- 
nomena from various fields of critical phenomena, namely magnetic phase 
transitions with dipolar anisotropy, biased percolation problems, and the mul- 
ticritical point of the Mott insulator to superfluid transition in Bose lattice 
models. 

Another variant of Amit and Goldschmidt's generalized minimal subtraction 
procedure has been used by Stephens and O'Connor [26,27] to investigate 
crossover phenomena ranging from finite size scaling problems to quantum- 
classical crossover. 

We close with emphasizing that the method could be applied to various other 
quite interesting crossover phenomena. Since it is quite likely that many exper- 
iments are not done in the asymptotic but in the crossover regime, it seems 
highly desirable to perform such calculations in order to get a quantitative 




(38) 



1/2. 
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understanding of the critical behavior. 
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